Phase Transitions in Confined Antiferromagnets 
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Confinement effects on the phase transitions in antiferromagnets are studied as a function of the 
surface coupling v and the surface field h for bcc(llO) films. Unusual topologies for the phase diagram 
are attained for particular combinations of v and h. It is shown that some of the characteristics 
of the finite-temperature behavior of the system are driven by its low-temperature properties and 
consequently can be explained in terms of a ground-state analysis. Cluster variation free energies 
are used for the investigation of the finite temperature behavior. 
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In recent years the theoretical interest on two- 
sublattice, uniaxial antiferromagnets has been renewed 
. This interest stems from the experimental work on 
Fe/Cr(211) multilayers by Fullerton and coworkers ||, 
where an antifcrromagnetic coupling between the Fe lay- 
ers is possible for a suitable choice of the Cr layer (11 
A). For an even number of layers, the spin-flop phase 
nucleates at the surface and, as the external field in- 
creases, this surface phase evolves into the bulk spin-flop 
phase Q. The small anisotropy-to-exchange interaction 
ratio that characterizes the Fe/Cr(211) system, makes 
this particular sort of multilayers amenable to the theo- 
retical modeling. Since this type of magnetic multilayers 
are isomorphic to the M11F2 class antiferromagnets with 
(100) surfaces, a classical one-dimensional XY model is 
adequate to model their magnetic properties at low tem- 
peratures . The early work of Mills Q] and Keffer ||] 
on the one-dimensional XY model, that lead to the iden- 
tification of the surface spin-flop transition, have been 
complemented and extended recently. The surface and 
finite-size effects on the ground-state properties of an 
XY chain have been investigated in terms of discom- 
mensuration transitions ^,|| and the analogy between 
the one-dimensional XY model and Frenkel-Kontorova- 
type chains has been elucidated M . These investigations 
have set our understanding of the rich magnetic behavior 
in Fe/Cr multilayers, where finite-size and surface effects 
are equally important, on solid physical grounds. It is 
appropriate to note that the passage from an inherently 
three-dimensional structure, such as the Fe/Cr multilay- 
ers, to a one-dimensional structure is based on the as- 
sumption that lateral fluctuations within each layer can 
be disregarded with respect to the interlayer fluctuations. 
The effect of confinement (surface plus finite size) in an- 
tiferromagnets for which the intralayer fluctuations are 
important has also been the subject of previous investi- 
gations |§-[ll|. 

In this paper we provide a brief survey of a recent study 
of ground-state properties of bcc films with [110] sur- 
face orientation fbl and relate the previously observed 



topological features of the phase diagram |^|,|rTJ to the 
zero temperature properties of the system. With these 
objectives in mind, we consider body-centered antiferro- 
magnetic Ising films with surfaces oriented in the [110] 
direction. In the (110) planes of a bcc structure, each site 
in one sublattice has nearest-neighbors in the other sub- 
lattice [lq ]. For nearest-neighbor pair interactions, the 
Hamiltonian is the following: 

H = Jb /J <7i<7j + J s 2_, aia 3 
ij^surf ij'Gsurf 



iebulk 
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where the spin variable cr, takes the value of +1 or — 1 
depending if the spin at site i is pointing up or down, re- 
spectively. We have assumed that surface sites, in layers 1 
and N for an TV-layer film, experience a surface magnetic 
field h in addition to the external field H . We can think 
h as the surface perturbation on a highly anisotropic an- 
tiferromagnet (Ising-like) slab, caused by the presence 
of ferromagnetic layers in a FM/AFM superlattice. We 
can also consider h as the wall-particle interaction in a 
fluid confined between two parallel plates, when the usual 
transformation pi = ^(l+crj) is used to cast Hamiltonian 
(|l|) into a lattice-gas model. The wall-particle interac- 
tion (h) is responsible for the condensation of the liquid 
phase at lower chemical potential than it is necessary in 
the bulk (capillary condensation) [^6|-^0|. 

Phase equilibrium in confined systems is very sensitive 
to the boundary (interface) conditions defined by the sur- 
face field h and by the surface coupling J s pi|^2[| . In 
this paper we specialize ourselves to the case of nearest- 
neighbor interactions and localized symmetric surface 
fields; that is, the field at each surface is the same and 
acts only at the surface sites [see Eq. ([[])]. In the re- 
maining of the paper, the effective pair interactions, the 
surface field (h), and the bulk external field (H) shall be 
expressed in terms of the bulk AF coupling (Jb >0). The 
ratio of surface to bulk coupling is restricted to positive 
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values and it is denoted by v. 

Even when h is zero (neutral boundary conditions), the 
disruption of the translation symmetry due to the sur- 
faces results in a "missing neighbors" field h m . The sur- 
face field h m is responsible of the inhomogeneities in the 
magnetization profile near the surfaces. When Eq. (0) is 
reinterpreted as a binary-alloy Hamiltonian, the missing- 
neighbors field, along with h, accounts for the surface 
segregation phenomenon, i .e., the enrichment of the sur- 
faces with one component [pO 22 23]. In the following we 



shall consider only the case of h > 0, since the results for 
h < can be obtained straightforwardly from the sym- 
metry properties of Hamiltonian (0), as it is discussed 
next. For zero surface field, Hamiltonian (0) is invariant 
under the transformations <7j — ► — <7j, H — > —H. For 
neutral boundary conditions (h = 0) the ground-state 
and the finite-temperature properties of the Hamiltonian 
are symmetric about H = 0. A positive value of h breaks 
this symmetry by favoring the spin- up states at surfaces. 
Both zero- and finite-temperature properties of Hamilto- 
nian in Eq. (0) become asymmetric with the applied field 
H. However, for nonzero h the Hamiltonian is still in- 
variant if we extend the above transformation to include 
h -» -h. 

The selective nature of the surface field has varied con- 
sequences on the properties of Hamiltonian (0), since 
the equilibrium states are defined by the competition be- 
tween the Zeeman and the ordering energies. The order- 
ing contribution to the Hamiltonian, first term in the rhs 
of Eq. (0), favors AFM structures whereas the Zeeman 
energy in third sum of the rhs of (0) promotes FM struc- 
tures. An additional Zeeman contribution arises from 
the surface field [last term in the rhs of Hamiltonian (jjj)], 
which competes with the surface ordering tendencies [sec- 
ond sum in the rhs of Q] to define the equilibrium state 
in the film. Thus, for positive and large values of H, 
applying a surface field is of little consequence since the 
stable state is one with spin- up at the surfaces. For low 
and negative values of H, where spin-down states are 
likely to occur, the surface field actually may give rise to 
an antiferromagnetic surface state. 

An analysis of the ground states of Hamiltonian p) sin- 
gles out the ground-state (GS) sequence in Fig. 0(a) as 
the stable sequence for large h [Q . The nomenclature in 
Fig. is as follows: the intra- and interlayer coordination 
numbers are represented by zq and z\ , respectively, with 
the bulk coordination number expressed as z = zq + 1z\ . 
The parameter z s — zqv + Z\ can be regarded as the sur- 
face coordination number but actually accounts for the 
surface energy [recall that all quantities in Eq. (0) are 
normalized to Label Tl / I / Tl stands for a iV-layer 
film with AFM surfaces and down magnetization in the 
remaining (iV — 2) layers. 

From Fig. 0(a) one can see that a GS structure with 
AFM surface coexists with a ferromagnetic bulk for 
H G (— z s — h, — z), while the contrary occurs for H £ 
(z s — h,z). In between, i.e. for H E (— z, z s — h), the 
GS is AFM in both the surfaces and the bulk. The film 
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FIG. 1. Schematic representation of the ground-state and 
finite-temperature phase diagrams for films under intense sur- 
face fields. As a function of the external field H , a film of iV 
layers transits between the ground states (GS) displayed in 

(a) if h < h g or in those showed in (c) if h > h g . For the 
former case, the corresponding H-T critical line is shown in 

(b) in a thick solid line, (c) For h > h g a ferromagnetic (dis- 
ordered) gap intervenes between Tl / I / Tl ground state and 
the zero-magnetization structure. The characteristic field be- 
tween the disordered gap and the GS with surface AFM order 
is H ag — z s — 2zi — h. For h slightly above h g , thermal exci- 
tations turns the ferromagnetic gap into a disordered region 
in the H-T plane [dashed line in (b)]. For very intense sur- 
face fields the otherwise connected AFM region splits into two 
separate critical curves [thin solid lines in (b)]. See the text 
for further explanations. 

displays an ordered, compact domain from H = —(z s +h) 
to H = z. In Fig. 0(b), in thick solid line, we show the 
critical curve (schematic) in the H-T plane associated 
with the GS sequence of Fig. 0(a). 

For negative values of the external field (H < —z), the 
surface field favors the AFM ordering at the surfaces but 
also promotes the decoupling of the surface layers from 
the rest. Therefore, the inner layers closely behave as a 
(N — 2)-layer film with neutral boundary conditions. In 
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Fig. |l|(b) with thin lines and appropriately shifted, we 
have plotted the corresponding critical curves for a 2D 
square lattice (left) and the corresponding (N — 2)-layer 
film at h = (right). Observe that the shoulder shows a 
maximum temperature ~ uT sur f, where T S urf is the Neel 
temperature of the 2D square lattice. 

The ground-state sequence in Fig. |](a) becomes unsta- 
ble upon an increment in the surface field, and the GS 
sequence of Fig. 0(c) is then adopted by the film. Note 
the a disordered gap (GS T / I / T) an d a new zero- 
magnetization ground structure (T/l/TI/l/t :an 
up-magnetization state at the surfaces, subsurface lay- 
ers down magnetization and the rest in the AFM state) 
appear in lieu of the homogeneous AFM-GS structure. 
It can be shown [[t4| that the transition between GS se- 
quences, from the sequence in Fig. Qa) to the one in Fig. 
||(c), occurs at a surface field value h g given by 

h g = z v + (z Q + zi). (2) 

For a surface field slightly above h g , the disordered gap 
transforms itself, via thermal excitations, into a disor- 
dered region in the plane H-T, right in the middle of 
the ordered region [Fig. 0(b), dashed line]. In other 
words, the phase diagram is composed by two critical 
lines [dashed and thick solid lines in Fig. 0(b)]. Upon 
high-temperature cooling and for H £ (z s — 1z\ — h, —z), 
the system undergoes a phase transition from the high- 
temperature disordered state to an AFM state. A fur- 
ther decrease in temperature drives the system into a 
low-temperature disordered state. 

Intense surface fields increase the disordered gap at 
zero temperature and, as a consequence, the height of 
the associated disordered region rises. At h — h s , the 
AFM domain splits into the surface and the bulk critical 
curves. In Fig. 0(b) with thin lines, the critical curves 
associated with the surfaces and the bulk are presented 
for h > h s . It is shown in Ref. ]IJ] that the splitting 
point corresponds to a saddle point in the Hessian of the 
free energy as a function of T and H [Q . We have used 
the pair approximation of the cluster-variation method 
(CVM) |2fj] to evaluated the finite-temperature proper- 
ties of Hamiltonian (|l|) . Previous work have shown that 
for nonfrustrated lattices, such as the bcc and simple 
cubic, the PA gives reliable results for the qualitative as- 
pects of the phase equilibrium in restricted geometries 

In a sense, the splitting value of the surface field, h s , 
represents at finite temperatures the role of h g . Both 
characteristic values of the surface field h g and h s , are the 
answer for the following question: How intense need the 
surface field be, in order to split the otherwise compact 
AFM domain, into separate surface and bulk ordered re- 
gions? At zero Kelvin, the answer is independent of the 
film thickness: When surface field reaches the value of 
h g = zqv + (zq + zi), the surface splits from the bulk 
independently of the number of layers. At finite temper- 
atures, the answer is more involved since now thermal 
excitations enhance the coupling between the bulk and 
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FIG. 2. Surface and bulk fields corresponding to the split- 
ting point h s and H s (left inset), respectively, as a function of 
the thickness of the film. A magnetization profile for N = 100 
film at h = h s is also shown in the right inset. Calculations 
were performed in the pair approximation of the CVM for 
bcc(110) films with v = 1. 

the surface layers. Figure shows two regimes of be- 
havior for h s as a function of the film thickness: for thin 
films the value of h s increases with N while the contrary 
occurs for thick films. This peculiar behavior of h s {N) 
results from the balance between the surface Zeeman en- 
ergy and the ordering energy. For very thin films the 
surface Zeeman energy easily overcomes the contribution 
of a bulk made of a few layers. In this regimen, increasing 
the thickness in the film is equivalent to enhancing the 
bulk contribution to the free energy. Thus, it is necessary 
to apply more intense surface fields to split the surfaces 
from the bulk. 

For very thick films the splitting value of the surface 
field shows virtually no change as the thickness in the 
film is reduced. Near the splitting point, the magnetiza- 
tion profile decays very fast towards the bulk state as we 
move from the surfaces to the inner layers (see inset in 
Fig. ||). The surfaces are too far away to affect each other 
and, therefore, h s corresponds to the semiinfinite value 
of the surface field hf. However, if the film thickness is 
reduced enough (N ~ 50 in Fig. |J), the perturbation in- 
troduced by the surfaces reaches the middle layers. The 
interplay between the surface and the finite-size effects is 
reflected as an increment in the value of h a as the thick- 
ness is decreased. 

In summary, we have shown that the rich magnetic 
behavior, previously reported in AFM thin films p,|TT|], 
is directly related to the ground state properties of the 
films. We focused on the thermodynamic behavior for in- 
tense surface fields, since in that case the otherwise com- 
pact antiferromagnetic regions splits into surface- and 
bulk-driven critical curves. In the bulk-driven critical 
curve, the surfaces are less ordered than the layers in 
the bulk. On the other hand, along the line of phase 
transitions driven by the surface, the bulk layers are less 
ordered than the surfaces. For surface fields such as 
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h > h s > h g , in which the critical curve is well sep- 
arated into the bulk- and surface-driven AFM regions, 
the surface order parameter vanishes with exponent j3 s , 
which in the mean field approximation used here, equals 
1 for the bulk-driven critical line, whereas [3 S = ^ for 
the surface- driven phase transitions. However, our pre- 
liminary results show that even for h < h g the surface 
exponent changes from j3 s — 1 to [3 S = i as the external 
field is varied from positive to negative values. The in- 
vestigation of the critical behavior will be considered in 
the future. 
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